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AB smog 

!She sul5critioal transonic flow past axisyraaotric and 
two-dtuonsional slender lo^es has heen siadied in the 
present work* She integral eciuation approach has heen 
adopted* For this nonlinear transonic small distur banco 
equation is transformed into an integral equation using 
G-reen*s theorem* !Ehe solution for the axial perturhation 
veloolty at a general point in the flow field is sought in 
the meridian plane hy diyiding the region of integration 
into rectangular elements, wherein the axial perturbation 
velocity is taken unifarm* fhus an iterative scheme is 
developed of which convergence is sought* 3?hG results are 
computed for a parabolio arc of revolution, parabolic 
airfoil and HiGiWOOdS airfoil* The results obtained are 
con^ared with experimental results and/or results obtained 
by advanoed relaxation methods. 
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Pressuro cooPfioient 

Goinplete elliptic integral of second Bind 
Gcmplete elliptic integral of second kind 
Bineness ratio of parabolic arc of revolution 
tunnel half height to body . chord ratio* 

laplaciaxL operator 
Mach number 
Porosity parameter 
Normalised oylindrical co-ordinates 
transformed cylinctoioal co-ordinates 
Normalised cartesian co-ordinates 
transformed oartesian oo- ordinates 
Ratio of specxfic heats ( s= 1 ,4 for air) 
Running oo-ordinatos in x,r,0 directions. 
Running coordinates in x, y directions. 
Perturbation velocity in x~diroction. 

Linear pc. turbation veloLxty 
thicloiess ratio 

Perturbation velocity potential 
iCransformed perturbation velocity potential 
Parameter satisfying Laplace equation 
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Xj JC j0 

x,S,9 
XjJf z 
Sfy,z 
5,P,V 

?,c 


Values of variable at different locations 
Partial derivative of variable v/»r.t. x,r,e 
Partial derivative of variable wl^*t, Xj^jS" 
partial derivative of variable w.r.t* x,y, z 
Partial derivative of varialiLe wa:.t. x,y,z 
Partial derivative of variable v/,r*t. 5,p,v 
Partial derivative of variable Vir.r.t.5,? 
Pree stream condition. 
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IJIST O P Flgg RES 

1* Slender body in axisyinmetric flow. 

2 . Shin a,irloil in tv/o-dimGnsional ilow. 

3 . Loxjation of grid point in r octangular jleBients . 

( axisyrimetr ic flow) . 

4. Looation of grid points in rootangular elements. 

(two dimensional flovO * 

5* Distribution of 0^ on parabolic arc of revolution 
(M^ t=f 0 *9, S'!! = iO, Dree air oaso) . 

6* Distribution of 0^ for parabolic arc of revolution. 

(M^ “0 *9, DR = lOj Dree air case) . 

7. Distribution of Op for parabolic aro of revolution, 

(M^ “0*9, DR K lO) wall interference case with 
porosity parameter 0*77 and wail radius 1,17, 

8, Distribution of Op on parabolic arc of revolution 
{ 11 ^ “ 0 *85, DR =s lo, wall radius i ,17) , 

9 . Distribution of on parabolic aro of revolution 

« 0.9, DR = lo, Solid wall at 1,17) , 

10 , Distribution of Op on parabolic aro of revolution 

,9, DR « weJ-l -- at 1,17). 

11, Distribution of on parabolic aro of revolution in 

£i?ee air at supercritical shook free mach numbers (DR = lo) • 

12. Distribution of on general parabolic- body of revolution 
in free air (M^ ^ 0 .9, DR a lo) . 
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15 , Distribution of Op on lUOA-oolS aicfoil at zero 
incidenoe in free air (M^ == o .V2) . 

14. Distribution of Op on 6 poroent tliiclc parabolic 
airfoil at z^o incidence (M^ ~ 0 »82c ) 

15 Distribution of 0^ on 6 poroent tld ole parabolic 
airfoil at aero incidence y/ith vyall interference* 

(M^ =s 0 .825, h/c S .0) . 

16. Distribution of Op on 6 percent thioK parabolic 

airfoil at soro incidenoe at tiacb number OtV and 0 »8 
in free air by present metliod. 





Despite imob advances in aerodynamic research for 
flight regime lying far beyond the so called sonic barrier] 
there still remains some probl^ois connected with the 
phenomenon that occur in the flight range straddling the 
speed of sound. Consequently this range transonic flow 
range - continues to engage substantial research effort 
and gained impetus due to possible development of minimum 
drag transonic cruise configuration for transport aircraft 

To simplify the study of transonic flow and for ease 
of building the mathematical model, following assumptions 
are made % 

i) the fluid is ideal and perfect gas 5 
ti) the medium is continuimi and homogeneous; 
iii) the flow is steady and Irrotational 
and iv ) the body considered is slender. 

Under these conditions with perturbation velocity 
^y5 being quite small the equation for perturbati 
velocity potential is given by Prandtl-Glauert equation j 


a • 111 ) 


0 
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The above equation governs subsonic and supersonic 

Tlow quite reasonably* However since in the transonic 

range 1 the coefficient of 0^ on L.H.S. becoraes very 

small hence an order of magnitude analysis of gas dynamics 

equation necessitates retention of some non-linear terms 

27 

also. The final equation takes the form s 

(1 - M® ) 

The above equation in cylindrical co-ordinates is t 

(1- ^rr ”2 ^G0 r '^r ” ^ i+i) 0 ^ 0 ^ 

r 

This small perturbation equation is equally valid 
for transonicj subsonic and supersonic flows. In the x-r 
plane, it is a mixed hypdrbolic parabolic and elliptic type 
nonlinear equation i.e. for : 


Cl - 14- 


(y4.i)/^3 


""V 0 elliptic (subcrltical flow) 

= 0 parabolic (critical flow) 

< 0 hyperbolic (supercritical flow) 


The shock free supercritical case is of particular interest 
as it gives minimum drag for the body. 


To solve the transonic flow equation many methods 
have been tried of which relaxation method is quits standard 
method but needs a lot of labour* An equally good technique 
is integral equation approach. It is less time consumiiig 
and attempts are being made to make it of similar accuracy 
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as tjiat of relEucatictti method* A brief survey of integral 
equation method is presented in the following section. 

1»2 Idteratune Sur^y i 

Integral equation approach is oie of the earliest 
successful method for solving the direct transonic flow 
problem. !Dhe main idea is to convert the transonic flow 
equation into an Integral equation using Green’s theorem. 

The solution to the integral equation is obtained by 
iterative technique after certain modifications have been 
incorporated into the equation so obtained. 

Among the most noteworthy earliest atteiai^ts for the 

solution of transonic flow problem by integral equation 

method are three related studies that deserve special 

mention. These include those carried out by Oswatitsch, 

Gullstrand and Spreiter and Alksne, and hence the approach 

1 

adopted has been termed as OGSA approach by Ferrari « 

2 

Oswatitsch deduced a ncn-linear integral equation over the 
whole plane flow at zero lift and made approximate calcu~ 
laticaas for some aerofoil in the lower transonio range. The 
work was later extended and iHustrated with applications by 
Gullstr^d ■ • more recent exposition of the method is 

Q 

given by Spreiter and Alksne , who worked out the method 
more thoroughly* Their results are fairly good for the 
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subcriijioal flow around circular arc airfoil at zero 
incid^oe * Heaslet and Spreiter^ applied the integral 
eq,uati< 3 iQ. njethod to transonic flow around slender wing and 
bodies of revolution* Special, attention was given to 
conditions resulting from shockwaves. Results are obtained 
for cme cylinders, wings and wi3ag body combinations at a 
free stream Mach no**«l and compared with experimental 
results • 

N^rstrad"^” extended the OGSA integral equation 

method to lifting and non-liftirg transonic flow problems, 

including the flow with shock discontinuities. She integral 

equation is transformed into a set of non-linear algebraic 

equations* (Do ensure a unique solution for supercritical 

shock free flows the method of parametric differentiation 

has been applied. His results indicate generally good 

agreement with experimental results and otherwise- accurate 
those 

results IjOce^/obtained by finite difference techniques. 

(Bhe integral equation method to solve the problem of 
high subsonic flow past a steady two-dimensional airfoil has 
also be^ used by Nixon and Hanoock^^, (Dhey have confuted 
the shook free two-dimensional flows around the lifting and 
non-lifting airfoils by solving the integral equations 
approximately, (Dhe results for two test cases show close 
agreen^nt with more exact numerical results of Sells^, 
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The method has also been extended hy Nixon to an airfoil 
oscillating at. low frequency in high subsonic flow. 

A feature common to all above integral equation methods 
is that the field integral is not evaluated accurately. In 
these methods using a functional relationship for velocity 
distribution, the double (surface) integral is reduced to 
a line integral in the plane of the airfoil. In the extended 
integral equation method of Nixon an alternative means of 
evaluating the field integral is developed. Ihe flow field 
is divided into a number of streamwise strips and the trans- 
verse variation of perturbation velocities across each of - 
these strips is approximated by an interpolation function 
in terms of values on the strip edges, Ihe field integral 
is then reduced to a line integral \'rhich is in turn evaluated '*• 
by quadrature. Ihe results obtained by using this method 
have been compared with those of the standard integral 
equation method (See Nixon ), It is shown that this method 
gives considerably improved results over the earlier integral 
equation methods. 

Recently Ogana^^ obtained results for two-dimensional 
parabolic arc airfoil evaluating the field integral quite 
accurately. He divided the region of integration into 
closed rectangular elements, and considered velocity to be 
constant in eaoh element 5 the integral equation is then reduced 
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to matrix eq_uation solvatle by suitable iterative scheme fo 
subcritioal flov^s # 

She influence of wall interferaneG effects on models 
tested in transonic wind timnels has been an open question 
for many years, fhere han been an increasirg need for the 
largest model possible* She relaxation method has been 
applied to the problem of wall interference effects. Bailey 
treated the slender bodies of revolution at transonic 
speeds in the presence of wind-tunnel wall boundaries, whili 

21 pp 

Murman and Kaopraynsloi ha've considered interference 

effects on two-dimensional airfoil, lifewum. and Klucker^® 

have extended relaxation, method to determire interference 

effects of transonic flow over finite lifting wings in 

rectangular tunnels. Murman, Bailey and Johnson^ have 

given an upto date coiiputer program to solve the transonic 

small disturbance equation for two dimensional flow for 

lifti]:g airfoil which takes free air case as well as variou. 

wind tunnel wall ccnditions. Ihey used finite difference 

formulation using iterative sucoossive lih^ over relaxation 

(ShOH) algorithm. ^ integral equation method has also 

25 

been used by Kraft for wind tunnel boundary interference 
at transonic speeds on thin airfoil in two-dimensional wii 
tunnels* However, field integral is reduced to a line 
integral over the airfoil by the use of arbitrary approxiiaal 
function and the results are obtained on an airfoil* 



Xhe integral equj'ition method applied to two^ 
dimeneionaX flow has provr^d to glvd good rc<=^ults, Howeirerf 
the attempts to solve ths problem of mixed tr?nsonlc flow 
past bodies of revolution by integral cion method have 
been dismal* 

In present work an iiitegi^al e.^via^tlon suitable for 
axioymmetrlc flow past slender bodies of revolution at 
transonic speeds is obtain^-d in the fre- air and in bho 
presence of porous wall boundarios* The numerical solution 
to integral equation is sought in the meridian plane by 
divxdjng the region of inbegrotion into rectangular 
elements wherein the velocity can be considered uniform* 

The pressuro field ib computed on and away from the slender 
body including force and aft body effects* to application 
to parabolic bodies of revolution with end vrithout sting is 
presented ab subsonic Mach number.'^ yieMing subcritical 
and slightly supercritical ‘^hock free flows. 

Similar analysis has oe^n carried ou^ for two-^ 
dimensional flows also , Parabolic airfoil has been 
considered in the free air and in the presence of porous 
wall* An attempt is also mads to determine prsbjure 
field for NACA OOlS airfoil at aero incidence. 
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PR-OBM aPBO Ig lOAO^IOM 
2^X Introduotxon i 

Most of the aerodynamic problems are related with 
finding pressure field over tv/o-dimensional or axisymme 
bodies j so tha'li one can predict coefficients of lift, dra^ 
or mcment over a body in certain situations* Here we are 
interested in finding pressure field on and away from body 

I 

She axisymmetrio problem is considered first and then the 
flow past two diaensional non lifting airfoil is treated* 

2 «g Axlsynanetrio Plow s 

Consider the invisoid, compressible flow about a 
slender, pointed body of revolution in cylindrioal co= 
ordinates x and r (made dimensionless v^ith respect to the 
length of the body) with x axis parallel to the body axis 
and free stream velocity (see fig* 1) • Assume she 
free flow so there is no vorticity involved* Under this 
assunption a velocity potential can be used to calculate 
the flow field. In particular a perturbation potential 
con be defined such that the perturbation velocities, made 
dimensionless with respect to the free stream velocity, ar 
h = the x«axis and v » parallel to r-ax 
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(The governing differential equation for transonic flow 
given ly slender body theory can be written 

(1- wf, )0xx + + ? I^r = o +1) 0^ , 0^ (2 .1) 

where is free stream mach number and y is the ratio 

of specific heats, 

Ihe term associated with 0 makes transonio equation 

A 

both nonlinear and of mixed e Hip tio-hyp orb olio type, Ihe 
mixed character of the flow field here may occur with 
local supersonic regicns embedded in a subsonic flow. 


lo complete the specification of problem* boundary 

ccnditicns must be given at the body and in the outer flow, 

Ihe flow tangency condition at the body surface given by 

the first-order slender body approximation is that near 
27 

the body axis j 


lim- r 0 
r “* 0 


R 


m 

dx 




( 2 , 2 ) 


where r = R-(x) defines body 

S(x) defines cross sectional area of the body and 


prime denotes derivative of function w.r ,t, its 
argument , 

I'or a body in free air the perturbation velocities 



10 


vanisil at in‘£inityj 

0 as ^ (2.3) 

^ X 

which is satisfied by setting 0 equals a constaiit, 
say aero, at infinity* 

in addition to the free air boundary condition wall 
boundary condition may be given to approximate inviscid 
flow in porous wall wdhid tunnel test sections, ihese 
boundary conditions can be used to illustrate wall-induced 
interference effects. Although the formulation is strictly 
valid for a circular test sootionf the results may be 
compared with a square test section of equal cross sectional 
area because the effects of wall interference at tie centre 
of tunnel should be relatively insensitive to the actual 
wall shape , 

Ihe average boundary condition for a porous v^-all, as 
derived by G-oodman^^ follows from Darcy's law for slow 
viscous flow through a porous medium. It is assumed that 
the average velocity normal to the wall is proportional to 
the pressure differenoe across the wall, which is a linearised 
approximation of Darcy's law for a thin wall, and pressure 
outside the wall is free stream pressure. v7ith the wall 
parallel to x«axis the porous wall boundary condition becomes; 

+P0^ — 0 at the wall r = rw (2.4:) 
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The q_uantity p is a porosity parameter or iv.oynold.s lumher 


of the porous raediuni, defined by : 


U, 


p = 


P k 


W t 


(g*5) 


Yifhere 

free stream "velooity 
free stream density 
k permeability of porous medium 

y viscosity of air 

t tunnel wall thickness 


k is detorrained by the structure of the porous 
medium and must be found e:xperimentally^^ • 

As porosity parameter vanishes i.e. p == o eq^uation 
(2,4) gives boundary condition for solid wall and the 
porosity parameter becomes very large i,e, p oo equation 
(2,4) approaches the open get boundary. Although truly 
porous walls are seldom used for transonic wind tunnels, 
the porous wall approximation is useful mean boundary 
condition for perforated walls, 

I’inally, the slender body approximation for the 

P’7 

pressure coefficient at points near the body is given by 

Op = - 80^ - (a.6) 

and on the body surface 

Op «= - -(<ffi./ax)^ 


(8.7) 
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1.3 !gwo-.Dijnensic3aal Plow * 

Here Oartesian ooordlnatos (made dimensionless 
w.r.t. the length of airfoil) has been chosen to describe 
the inviscid, compressible , shoch free flo?/ over a thin 
symmetric airfoil. Ihe x-axis has been taken to be parallel 
to airfoil plane and the free stream velocity (see fig .2). 
The perturbation velocities made dimensionless w.r.t, the 
free stream velocity are u = parallel to the x-axis 

and V ™ 0 parallel to the y-axis , 

V 

TIk governing equation for transonic flov/ 

27 

given by thin airfoil theory can be written as 

( 2 - 8 ) 


where symbols have same meaning as in previous sectiont 
The boundary condition for thin airfoil is 


given by 


27 


lim- 

y 


0y(X, +0) = 


dX + 


where y = 


+ 0 dx 

defines upper airfoil profile 
^ ’f-(x) defines lower airfoil profile 


(2.9) 


For symmetric airfoil t+(x) - «Y— (x) . 

For the airfoil in free air perturbation 
velocity vanish at infinity which is satisfied by setting 

0^ , 0^ 0 as y^ (2,l0) 


OQ 



5*or woU antcrf rcncc. cace the effect 

porous wall, parallel to free stream oan. he taKon by 

26 

applying porous wall boundary condition given by 


+ p 0^ = 0 at y ^ + y\Y 


where + sign is lor upper wall, - sign is for lower wall, 
and p is defined by (2,5). 


Ihe thin airfoil approximation for pressure 
coefficient for the point near the airfoil is given by 

Cp = 2/lrM^{[l +1^ m2 (1.(^4,^) 2 _ 2,] 


( 2 . 12 ) 
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UTWRAIi 


RffiJTHOD OR SOLUTION 


5 ,1 Introducti-on t 

In this chapter the governing nonlinear 
partial differential eq.uation of transonic flow is trans- 
formed into an integral eq^uatlon using Grroen’s Theorem for 
subcritical and shock free supercritical speeds (M^ < 1) » 
Axisyimetric and non-lifting tv/o dimensional cases have 
been dealt with* The axisymetric case will be treated 
first * 


letric flows 


To apply the Greenes Theorem to the axisymmetric 
transonic flow equation (2.1) one stretobss the coordinate 
system suitably. Here we make the following transformations; 


X = X 


r « pr ; 0 


where 


P = ^ 1 - 
K = CY + 1] 



(3.1) 


Now the governing transonic flow equation for axisymnetric 
flow (eq« 2*1) can be expressed as 
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^5aE + ^rr+ f C3.2) 

SimilaJTly iDoundary conditions eciuations (2,2), (2*3) and 
(2*4) beocme 

<’•=> 

0 at infinity (3 ,4) 


and for wall interference case at the wall r = : 

0 ^ + P/P ° 

How we shall apply G-reen*s Ihsorem to the ecLuation (3.2), 

Ihe theorem states that the following relation holds between 
any two arbitrary functions f and g having continuous first 
eind second derivative in the volume V bounded by the surface 
S 


IlfZsHt) - fi(g);]aT = //(f H - s H ) as (3.6) 


where n is inward normal drawn to the surface S and L is 
Laplacian operator in cylindrical coordinates i.e, t- 


L = 



1 

P 





here ift^) s,re running coordinates of domain T and S* 
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!So obtain, the solution ion 0 one ohosess 


f (? jrPfv)” ^(^»P>v) 


=* 0 (5 > p ) beoause of axisyranetric 

nature of 0 


s(5 » P :> V ) ™ ip ( e » P » V ) 

where i|» satisfies ) “0 which gives the fundamental 

solutions 





v/here H s= [I (x « C )^+ r^+ 


Z rp 


cos(e-.v ) 


Vg 


(3.7) 


fhen Equation (3 .6) gives the integral equation s 

H - I- (3 .8) 


Since ip and 0 do not have continuous first 
and second derivatives everywhere in the flow field hence 
we shall iUTrestigate to v/hat domain above relation holds. 
Ihe volume is (unhatohed portion in fig. 1) infinite volume 
Vgj (enclosed by wall if wall is present) excluding (i) a 
spherical point cavity at (5“X,p~r,v-e) where ^ 
is singular and (ii) the volume v/hich the body- 

covers. fhs surface S bounding the volume Y consists of 
i) surface S^i at <» for free air case or at wall 
if wall is present ; 
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ii) surface Sj, of spherical point cavity at 
(S==x,p^ r, v = e) 
and iii) surface bounded by the body. 

Now we cons id jr the volume a.id surface 
integral separately, 

5.S.1 Yolunie Integral s 

Nrom equation (3*2) 

K^) = \ 

1 
Wtm 

2 

where u = ^ 

hence the volume integral in equation (3,8) be comes a 

ly = I /// f dV (3.9) 

here Y « Vp - 

IDhe cavity volume Vp can be included in Y 
as it gives negligible contribution in the limit to (3,9) 
as shown below » 

Let us assume the radius of spherical cavity 
to be G I e “* 0 , we get 

4 3 

dY « Yp - ^ tc G 




is 




# 

• ♦ 


= finite » G^(say) 
X 


ly - g Iff ih dV 




X 1 ^ 4 „3 

2* ^1 S' ® ® 

f a,e^ 


- 0 as G -* 0 


hence we get 

V = Vj, - Yb 

How integrating equation (3.9) hy parts in g direction 
with houndary condition (3 .4) of perturbation -velocity 
vanishing at g - + «• one gets 

ly « " I (5.10) 

5 .2 *2 Surface integral 4 

# 

i) Integral over i 

For free air case the integral over the 
cylindrical surface at infinity can he eij^jressed as 

\ in - I- in 

« Sn: 

«» r/ f - ill 0 ) p dv d^ 

•noo 0 p b 

lim p -* 
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^This integral, wither given by equation (3.7) and with 
the assuiupticn that ^ ^ for e > 0 at infinity, 

vanishes • ^ ^ 

For wall interference case ^he integral over 
the finite cylindrioal wall surface 3 can be expressed as 


Irt — Iq *5 *“ 

^3! \ B. 


IS (^ 4 .^ 


w 


where r^ is the radius of cylindrical wall 


or la = 

W 


we get 

^w 


^ p ) dv 


o 


p=r. 


V/ 


Applying porous wall boundary condition (3.5) 




„ dv d^ (3.11) 


'w 


After simplification equation (3,11) leads to 
following line integral (see Appendix Ar-I) 




f + B) ?( (5 , sp - 0 I ^(5 , r„) 3d5CS.12) 


w — «> 


? E- 


where 


A 


- s B(b/ 2 , k) 


'W 


% 


(a.-b) "fo- 


K _ 1. Hli/Sf k) - B(V2, k) 

^ fa ~ 
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0 = Ee 

a =: )^ + (^ + ■>3 = 4 r r^, 

k = f B7a 


(3,13) 


E and i* being the complete elliptic integral of second 
and first kind respectively. 

Hence surface integral over can be expressed as 



where 


« 6 

6 


/ [[[.(A+B) 5 ) r^) ■* G o ^ ( 5 , r^) 3^5 

,oo r ^ 

(3 ,14) 


0 for free air case 

1 for wall interference case 


ii) Integral over Spj 

She point cavity is the sphere of radius G» 
G «> around the point ( 5 = x , p = r, v = e) the 

integral over Sp can be written as 



= C I - ^ S 

« illa|4 « -A— M 3 4m: 

^ 4 icR^ W aR 

s S(x, r) - G ^ (since here I^g) 

eR 

sa -* .? (x, r) G ** 0 


(3,15) 
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iii) Integral over Sg l 


For sufficiently slender body we approximate s 


lim 

P “*0 


- L. 

dn ■* ep 


hence we geiis 




p dy d 5 


(5,16) 


where £ is the non dimensional length of the body* 

^ p y 

In the limit p o, ^ is of the form 


iff S’(^)^T^P+g(^) 
Also frcm (5,7) % 


lim 

P -* 0 



r cos(0-v ) 


(S .17) 


Hence in the limit as ^*^0, and the integral 

(3,16) reduces to 

^ 2ff 

I«“ «/ / »{f|^pdvdC 

00 ^P 
lim p ■♦o 

Applying the tangency boundary condition (5*3): 

T . _ /f K 

0 -0 2 % 
lim P o 



2S 


JL^ f lO dc 


(3 ,18) 


Kow using (3,1j 0), (3.14), (3.15) and (3*18) in G-reen's 
theorem (3.Q) we finally obtain, a norw-llnear integral 
eq. nation for 


0(x, 5) 


_ K /._§a£j — 

° ^(£c )^+ 5*^ 

oo 

+ 6 / C(-A-+B) ^ (5 , r^) - 



G I 


+ I /// 't'r (C ,P )d-V (3,19) 

^ V ^ 5 


How to obtain an eq.uation for the axial perturbation 

Telooity we differentiate equation (3,20) w.r.t. x 
(see Appendix Ar4LI) and we get s 


u(x, r) w 



a 


/ 


S»fC)(S-£) 

CCS; -e )t 52 3V2 


d5 


+ 6 / IIa + B + u ( ? , r^)d^ ■ 

-,»s 

+ i //n s U-® dv (S.SO) 

P X C 

where 3) = r A 

^ r^ r 


and A, B' stand for sauie notations as defined in (3,13) 


4 



!Ehe equation (3*E0) can be rewritten as 

•* M oo 

u(5,r) = Uj.(S, r) + (A+BfD) u(C ,r^) d? 

1 

2 '^ 1 / X ^ > p ) 


V 


(3 .21) 


where Uj^(x, r) « ^2 Uj^(x, r) 

P 

n - 1 r s‘(e ) (x-e ) . r 

Uj^ being the linear value of perturbation velocity 


3 .5 !Bwo~ Dimensional Plow 

In previous section we derived axial velocity 
perturbation for axissnumetric flow. In the same manner 
the axial velocity perturbation in two dimensional flow 
over thin airfoil at zero incidence has been derived in 
this section. 


We first introduce the similar transf ormaticn 


• 0 • 


xaxj y~py?^ 



(3.22) 


where p s= "^1 « K = (Y + 1) Klf, i 

GO 4MW. 

Using (3,22) in the governing transonic flow equation (2.8) 
and boundary conditions (2.9), (2.1o) and (2.11), we obtain 



% 0-) - 


where X, (3c) ss 
+ 




\U) 


(3£4) 


0 at infinity 


(3.S5) 


and for the wall interferenee case, at the wall y ^ + y 


% = + I % 


(3.26) 


For tv/o«diinensional case the Green *s theorem which relates 
the surface integral to line integral is : 




where f and g are arbitrary function having continuous 
first and second derivative in surface S - S( 5 , C ) 
bounded by curve 0 = 0( C , ? ) and n is the inward drawn 
normal to the boundary. L is the Laplacian operator given 


,2 - g 

& , 9 

a ? a 


To investigate the solution for 0 , we choose t 

f •= ^ (e , o 


and g 
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where satisfies the Laplace equation L(\p ) = 0. !lhe 
fundamental solution of this equation is J 

c 5 In R (3,28) 

where R = h? - S)^+ (? - 5)^ 

Henoe equation (3,2V) together with equation (3.23) gives uss 

//* I dO (3.89) 

S '’■’0 

Nov/ we shall investigate the surface S and hounding curveC 
wl^re Iff and ^ have oontinuous first and second derivative 
keeping ir. view the houndary condition, 

such surface is infinite plane (enclosed hy wall 

if wall is present) excluding (i) a oircular point cavity 
Sp at (C == X ,C = y) an<i (ii) airfoil surface Sg • 

(see fig * 2) • 

$he curve 0 consists of (i) curve 0 ^ around infinite 
plane (or wall if wall is present) (ii) Ourve Op around 
circular cavity and (ixi) curve Og around airfoil surface , 

(see fig , 2) , 

Now we shall analyse the surface and line integrals 
of equation (3 ,29) • 



5 «3 Surface Integral; 


given iiy 


Frcm equaticai (3 ,29) integral over surface is 

Is = //'" 4S 

o 

= i " 0 ^ ^ ^ *30 ) 


}aere S = S^, - Sp - 


She cavity surface Sp can be included in S 
as it gives zero contribution in the limit to (3,3o) as 
shov/n below# 

Let radius of circular cavity be e, G ** 0, 


we get-5 

dS = 4u: 

u^ =! finite = Cg(say) 

= I (In e) 


(since here R = G) 


=5 0 as G 0 - 


hence we gets 

S = Sg, - Sp 

How integrating ec^uation (3,30) by parts in 5 direction 
with the boundary condition of perturbation velocity 
vanishing at infinity we get s 



2 ? 


5 ^5 Ijjie ^te^ral j 
i) Integral over 0^ : 


Slnco all perturbation vanish at infinity henoe 


for free air case integral over 0^ gives izero contribution 
while for wall interference case integral over 0^ gives a 


finite value 


In given by 


(see App endive B) 


I 



* _{ ( « . y*) { ^ C * C ? =■ yp + 'I' ( c =-y„) 3 


- / ^ = yw) ” ~ (3^2) 

0 for free air 

where $ = 

1 for wall interference case 
at upper wall 

and C ^ "" 

— y^ at lower wall 


ii) Integral over : 

Integral over curve Op is given as follows 


I 



n 


- 'P 0ji) 


« 110(5, y) 




1 2 tc 6 


where g o is radius of circular cavity Gp . 
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again In » 2 tc ^ (x, y) « G In G 

^ aR 


» 2ic ^(x, y) 


(since here Hr: g) 
as G ^ 0 


(3 *33) 


iii) 


Integral owr Og : 

Integral over "body curve Og is given by t 

dO 

■ n ■ '’ll' 


\ =/ i2f„) 


fX ^ Cx - Kf + (y -C A?, 


1/2 


0 


« A 0 


(G + ( C Jy)^ ^ ^0 


}:^re ^ ~ 0 ^ (x, 0"^) ~ 0^ (x, O”) 

= 0 (x, O"^) - 0(x, 0“) 

= 0 for synitTietric airfoil at aero incidence 


A ? 


from boundary ccaadition (3,21) we got; 

A = 


dl^c 5 ) a?_( 5 ) 

._ _ — 


+ 

T, 


= 2 


dx^ 

ST" 


for symmetric airfoil at sero 
incidence , 


which leads to 

1 . 5 1/2 dy.(5) 

- 2 / In []] (^ -x) + y H d? (3,3 4) 


B 


dK 



GomTjiniJig (3,32), (3,33) and (3,34) in eq.ua'tion 

(3,29) we get integral eq_uation for the transformed 
velooity potential as 

-1 i p 9 dY^(? ) 

y)= S) + y^H ^ ^ 

l«4aOO Q 

- f Di; (t = yw) “ * 5(5 =- yw^ 3 ‘15} 'I 5 

- ’*'5 “ 1 ^ '15 


Do obtain integral eq, nation for axial coB^onent of pertur- 
bation velocity we differentiate equation (3,35) w,r,t, x. 


5hus 


ti(S, y) 


CI5-e)®+'^ 


<1^ ( 5 J 1 

^ K 


dc 


§™ 

2)5 


CO 



^ ( c . yw^ { 



( c =yw^ 


+ 'I' (« = - Jw) 3 -I'i; ( 5 = yw> +'l'r^ 5 =" 
- {-fh S ^ 

o ^ 

00 

or u(z, y) » y) " 1^ / ^ 

00 P ^ 


(^==-yw^l 

^ 

4)t S 53 : 


(3.36) 
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where (x, y) = y) 

a^d VLjH, f) ae (3*5V 

^ ci^To^n 

here Uj^ is linear value of perturhaiion velocity for 
symetric airfoil at aero inoidence • 
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TOIERIO AL SOI EMB 


4 ,1 I: troductlon i 

In this chapter a numerical scheme is presented 
which has been used in calculating integrals involved in 
coGputation of perturbatiega velocity* Ihe integrals have 
been suitably expressed in summation form approximating 
velocity perturbation to be constant in small intervals. 

In this manner a set of algebraic equations are formed 
which are solved by method of iteration using suitable 
convergence criteria* 


4 »S . Axisymmetr 


ic Plow 


i 


We consider equation (3 ,21) for numerical 
evaluation. general axisymmetr ic body can be 

calculated numerically as follows . 


% (S, r) 


a 



_S1Ll JU i: 


'5/^ 


d 5 


which on integration by parts give 




V( jl-.x)^ + r^ 


I 


f ae 

V - - 




It- 


I ~ 






S^(x) In 



4x(£ «x) 




On the body surface r = E-(x) or r = pR eciuation (4,1) 
for pointed nosed bodies simplifies to 


u^(S, 0) 


1 



I t 


S (x) 


In 


ri ^ 

4x( i - x) 



(4.2) 


Ihe integral in eq.. (4.2) vanishes for ^ = x, She simpson's 
rule can be used to evaluate this integral quite satis- 
factorily • 


For the evaluation of linear perturbation 
velocity away the body, following procedure has been 
adopted* 


Ua 


(x , ?) = 


i- 2; s (Z- ) / 

^ -.- 1 . C(S-0‘+ 




- 5 


i:=l 
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(4^) 


where Xi - - 0 

fibre the body length has been divided into numbers of 
intervals of width 



Before aonsidering ectuation (3.22) wo attacli a 
] to X and P to r to sxniplify the analysis further. 
second expression on RHS of equation (3.21) can then be 
expressed as « 


CO 

6 / (1 + B + B) u ( ^ , r^) d^ 
n 


^ 6 


i==l 


'’1vjS5 



(4.4) 


where E^) 


n -^1 


'’iwSE 


f (A + B + B) a; 




Here the wall has been dlTidcd into n number of steps of 
size 26.: ar.d the Yelocity at the mid point of each step 

pi- 

hag been assumed censtant over that step , 


3Jo sirtplify further we proceed as in ippondix G 

and finally obtain t 




w 

l+lc(X^)’ 


(4.5) 
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where = 

Xg = 

^ = 

lo evaluate the volume Integral of equation (3.21) 
the region of integration has boen divided in reotangular 
grid in meridian plane , Ihe computational domain is defined 
by gi?id network as shown in fig. 3. Ihe perturbation 
velocity at the centre of a grid has been taken constant 
throughout that grid. Hence we get the volume integral 
after sin^lification as (see Appendix D) 


n + q 


X 


S 


i i 


Sn - X. 
o iX 








I /// f - ^ (( ,p )a.T 


n m p 
i=l 3^1 


^ijSP 


(4.6) 


where u. . 

and 


u(Xj_> 


r .+ h. 

^ / ' ■ / 


/ 


2« ^i+ h 


Si 


^ „ p d? 

X 


dv d p 


“ ^(2-]^)) 3}(Z]^) 

- &(l2) E(ii:/ 2, k(Xg)) $(Zg)3 (4*7) 


,** 


^ * 
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where G(X) 


^ r+(rp+ rp2 


1E(X) = 


k(X) = 


tan“^ ^.i” - tan”^ ^,1“ *~ ^n 


.1 - tan 

X 


X 


4 Tp r . 

^^■ipii^i^n^ >ijln ji 


T ?r^ir7ix 


Pinally using (4*4) and (4,6), etiuation (3,21) can he 
ej^ressed in iterative form as 


«(^) ^ 

^ + 6 y, 1^1 


_^(n-.l) n m 2 

“iwSP "iw ^Ijsi? hi 


(n-1) 


(4,8) 


g (say) 


Y/here 

li 

d (Xg, Tp) 

and 

U 

•M / »* V 


SP 


Pc..' convergence of abjve itcrati-ve ec^nation ^ u, 


has been added to both sides where x 
so that we obtain 
(n) n 


SP 


3U 


'SP 


(n—1) n m 

'SP i=l ^ i==l 3=1 

.(n-1) 


^P ” C "^iwSP '"iw 




^ ^rfn 3 / C 1- 2a„^.^ u 


_ (n-1) 


SPSP SP 


"SPSP '"SP 


n 


(4.9) 
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3!he convergence is achieved provided a^pgp ^ O.S/Ugp 

d *3 !Ey;o -. dimensional Plow % 

32he numerical iterative form of the integral equation 
(3*36) for two-din^nsional case can be obtained in the 
similar way as in § 4 ,2 , 


l^or most of the two-dimensional airfoils analytical 
results are available for linear value of velocity potential* 
fhese results for parabolic airfoils and ITAGA OOXI airfoils 
are presented in Appendix E-II . 

fhe wall integral in equation (3,37) is numerically 
e:^ressed as 

CO 

“ ^ -X 3 (5 , yp { ^ C'1’5 (C = +1'5 (c =-%) H 


- (c = y„) + ifjCc = -yp) <3? 


n 

i=i 


^iwSP 


u, 


rw 


where « u (x^, y^) 

introducing x - Xg and y - 



"’1'^ (c " y^) 



3*7 


~ If ^ f = yw^ + n 




X.- 6^ 


x.+6^ 


i— *-1 C 1 3 Cq-^ 5 

+ L^!“i _-^ + tan- J—— ;] } 

^v/"*" %> x.«6. 


yw - yt 


1 1 


1 ,p 


" Si~^2p' C 


in 


( yw"* yp ^ "** ^ ^ yw"^yp ^ ^ 

(V yp)®+ ^ " (vyp>^ + 4 


+ 


-.1 - - ^ 


yw " yp 


- tan 


^ ^1 . a.„„-l ^ 


+ tan' 


yw- yp 


yw+ yp 


1 ^1 

tan--^ ~=^i: ■} 

V yp 


where = Xg - + 5^ 


""S - 


Before representing surface integral of equation (3«36) in 
suitable form for numerical evaluation we use the symmetric 
property of 1i( 5 > ? ) for symmetric flow and limit the 
integration region to upper half plane only, as follows? 


" ^ 


1 


/ f 


y^ (jP 


' Xq ) "^yp) 


-Jrp C (§-Xg)^+ (^-yp)^Il 


-.2 




J. 

4n: 


^ p « p (ym is extent of C) 

^ ^ t rr. ! T, \(^ -L 


j” / r -(c-yp) 


•oo o 

- .2 


C(5-%)''+i:?-yp^‘^ 

Cc— Xg) — Cc ^ yp)^ 

c7e-5s)''+ fetyp)‘=3^ 


+ 


c[(e,Od5 d5 
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2 (say) 

Now the region of integration is di-vided into rectangular 
elements, as shown in Pig. 4, She perturbation velocity 
at the centre of each element is assumed constant within 
that element, Ihis gives s 


n m 


2 


13 


^IjSP 

_2 




a. 


IjSE' 


X4+ >4 .0 o 

1 r ^r(5-Xq) - fe-yp)® 

/ / 


V -*-0 • X. — 6 > ^ 
*^3 ^3 ^ i 




(e-xg)^ - (e-f y-^f 
[I ( 5 -Xq) + (c+ yp) I] 


d5 dc 




^ f 


yj+h. 


Xg « c 


Xg - c 


C|'-Sg)^+ C’S-yp)^ (C-Xg )"+(<?* yp)''J 

% 


-“TI2 


V ^i 

dc 


^3^ -3 - 




^3“ 


PM . 

^i*^i 


X. 


;r;F 


L ^ + (? +C5+yp) 

s, Zp 




dC 




-1 4.._-l V 


15. ^ j *b8xi 

y-D- c -~j ^3 ^3 


X-, 


^ .,-1 ^p“ ^ 

+ tan 

r. 


^3 ‘^0 
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ra — Q, ^ + tan**^ -i 

% X, X, X, 

-1 "^4 -1 1 x„ . T, 

«• tan Y“ ^ 4 " tan"' - tan**^ “1 

^1 ^1 Xi ^ ^ 

= yp - + aj j Xg « yj + y . - q. 

Tg sa y^ « y^ *- j ^4 “ Xp + Y ^ + 

2hus equation (3 .36) in iterative form 1300011163 

^Cn) ^(n-1) n ^ (iu.1) n m g(3^^) 

^SP “ ^SP * hwSS \w ^ijSP ^ij 

( 4 .I 0 ) 

where ^ tt (xg, yp) 

“igp = ' %'> 



SiMEJPHLS, 

MD, piSCTJSSfflH 

P rvL. li^'troduQtiQn j - 

In tflils chapijex 3resu,li:s of QCtupU'ta'fciuns arc presenijed 
and discuss ed* Resullis ai?e compu'tred at zefo angle of afiaclc 
for parabolic arc of revolution, parabolic airfoU and 
I^AGA^OOlS airfoil* Ihey are compared with ciKperimental 
results and/or that obtained by relaxation tochnipue* Ihe 
discussion for axisymmetric flow and two dimensional flow 
are dealt with separately, listing of the computer programme 
is given in Appendix F, Ihe oomputations were carried out on 
DEG«1 o 90 System. 

5,2 AKisymmetric Plow s- 

Yarious results for parabolic arc of revolution of 
fineness ratio lo v/ith sting at 0.854 have been obtained at 
different Mach numbers in free air and in the presence of 
wind tunnol wall. Ihe aft and fore effects are also included. 

For free air case integration was oarriod out in the 
range-0 *28 < x < 1.16 and 0 ^ r <1,1, Ihe stroamwise grid 
sise of A X = .04 was chosen, She transverse grid size A„r 
in physioal plane was ,2 except on the body surface v/here it 
was rl» Overall 37 x 6 rectangular elements were taJeen. 

Fig. 5 shows distribution of Op on the body at Mach 
numbac o *9 . It agrees quite satisfactorily with Bailey^^ and 
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experimental results of laylor and McDevitt^^ except in mid 
region* iBhe results for this case are also presented away 

from the body and conpared ?/ith experimental results^® (see 

fig. 6). lEhe convergence is achieved in 4 iterations'. 

Pig* 7 shov/s a comparison between present result with 

porous wall at radius 1.17 and porosity par^mieter 0.77, with 

29 

the experimental result ^ carried out in similar conditions^O 
Ihey are in close agreement exccept near the sting. Ihe 
integration was carried out in the range 0< x < 0.854 and 
0 1 3x < 1.17 without taking fore and aft effects at Mach 
number 0.9. It took 14 iterations to converge which indicates 
the porous wall interaction v/ith flow perturbations. Ihe 
results are plotted on and away from the body. 

Pig. 8 shows a comparison between three oases i) body 
in free air ii) body in presence of porous wall with porosity 
parameter 0*77 and iii) body in presence of solid wall, at 
Mach number 0.85. She axial integration range was - .28 < x 
< 1*16 and the wall v/as taken at radius 1.17. Par conver- 
gence the number of iterations were 3,9 and 5 respectively. 
2he porous v/ePl case do not have much difference over free air 
case. In the solid wall case the flow shows increased accele- 
ratim in the mid region. 

A comparison betv/een free air case and solid wall case 
ail- Mach no *.9 w'ith fore and aft effects is shown in fig* 9. 
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More iBcr eased acceleration oX the llov.r is observed here 1 
mid region v^hioh signifies that with increase of Mach 
niimhar the effect of wall interferenoe increases. !Bhe 

convergence v/as achieved in 4 and lo iterations respective 
for the t?/o oases , Porous wall interference calculations 

did not converge in the limit of 15 iterations. Hov/ever, 
without for and aft effects it converged in 14 iterations 
as mentioned above. Oorrespondingly solid wall and free 
air calculations were done without fore and aft effects 
(see fig* lo) * Ehey took 8 and 4 iterations respocti-yoly » 

fhe prestmt method wortn also for shock-free super*- 
critical case at Mach numbers 0.950 and 0.955 for free 
air with aft and fore body effects . Gonvorgence v/as 
achieved in 7 and lO iterations respectively. The pressin 
distribution on body surface is presentod in. fig. 11. 3}he 
result of Mach mmbai: 0.950 shows agreeable matching with 
esjp erimental results^^. 

Pig. 12 shows a saniple calculation for general 
parabolic body of revolution having maximum thidoness at 
50 perccait and 70 percent of body length. Oonvergcnoo was 
achieved in 4 and 5 iterations respectively at Mach numbers 
0.9 with aft and fore effects. 

3!he QPIT time varies from 40 secs to 1*2 minutes doper 

on number of iterations taken. 
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5 #3 ^ ;gviro..]>imeiiB3.prtaI ^lo v/ 

She ccaaputatiiorLs were carried out for NAG/w-OOlS airfofX 
and. 6 perceni; 'thick parabolic arc airfoil* 

pressor 0 coefficient over NAOA-OOlf airfoil was 
obtained at Mach number 0 in free air * fho integration 
was carried out in the range -.28 < x < 1*28 and *12 < y 
<9,12 and 29 x 6 rectangular elements were chosen with 
streainwise grid sise A x “ .04 convergence v/as achieved 
in 6 iterations talcing 2 minutes of GPU time, She results 
Y/ere cempared with lock's finite difference numerical 
results . Phey match well nesr the peak, hov/over, present 
is slightly lower negative at other points (see fig. 12). 

Results of 6 percent thick airfoil is given in fig .14:, 

15 and 16. She integration v/as carried out in the range 0 < x , 
< 1 and -4.56 < y < 1*56. 25 x 6 rectangular elements 

were chosen with stresmwise grid siae A x ;= .04, Rig 14 

IQ 

shows a ooinparisoo hetween the result ottained by Steger 
and present method at Maoh no. .826 . Result for present 
method was found slightly lower negatiTC than that of Steger' s. 
Oonvergenoe was aohleyed in 5 iterations taking do sees of 
OEU time. Sis- 15 shows the effect of solid wall on tne 
pressure ooeffiolent on the body at Maoh number .825. Ihe 
wall was taken at (h/O) = 2 .0 and range of totegration was 
0<x<l and_2<y<2 with 25 x 6 grid points. Ihe 
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solid wall interference gives higher negative pressure 
coefficient in the mid region of the body with percent 
increase at the mid point, Ihe oonvorgence v/as achieved in 
9 iterations Y;ith GPU time of 72 secs* Fig. 16 shows 
pressure dd^tribution in free air for Mach numbers .7 and .8. 
She converged in 3 and 4 iterations respectively taking 
23 secs and 31 secs of GPU time » 

5*4 Qonclusion 

An integral equation approach has been described for 
nuiaerloally cal-oulating inviscid subcritical -transonic flov/ 
about slender bodies. It has also been demonstrated that 
the present technitde can handle shock-iree supercnrxtical 
case. Ihe accuracy of the results obtained by the present 

technique is almost of the same order as those obtained by 
relaxation procedures . 

Oaloulations with wall boundary conditions have shown 
applicability of integral eq,uation toohniciue to the study 
of wind tunnel wall interference , 

Further this technique can be applied to lifting 

problem. 
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A-I Simplification of expressicin. 


E^^ressicaa (3 .11) is an integral over wall gi^ven by 

CO 3^ 

[o = - / / ('f' I 


4- u Sn; - ” „ 2n: 

f 0 f 4*1 «r / ^ /4 *q1'^15 

“ —CO ^ 0 0 ^ 


r p « 
w 


s^ijf p 


'W 


(a«i) 


1 

H0v7 consider each integral over v .sepeirately with 4!== 

1/2 

H = [^(x -c)^+ ^ 2r r^ oos(e- v)^ 

Without loss of generality we tato e « 0* Hence we obtain 


2% 

f 4j d V - 


1 f 11 


0 


L. 

% 


R 



(A^) 


^(x-O^ (r„, + 1 *)' 


/4B 

where 3c = / “-X/ 

\ (S-?) 3 ^) 

"see 291*00 of refer-noe ^ 


Again, 


1 

2-ji 

/ 

0 

f - r oosv 
w 

dv 


/ tKp<l'>=- S 

0 

""F 


It 

1 

C(^W 

1 

HI 

% 

*m r 

+ r j 

0 

dvQ 
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(see 291i0l and 291*06 of reference 


Ifenoe e^^ression (A.l) leads to ; 

B( it/S, k) 








. F(n/2, k) - E(7i/E,1c) , - , ^ , 

* k .r ' . n-^ ^ >V 


f 5 


E(n/2,lc) 




(S-5)^+ (", + 


I'X" 


■] 


de 


(A *4)) 


A-II Bifferentiatlon. of equation (5 *19)... ^ 

Dbe differentiation of equation (3.19) w.r.t, x is 
straight forward except the wall integral term 

!<. = / CC^+s) V- °fh 

% -.09 

On differentiating (A. 5) w*r,t* x , we have 


(A*5) 


C90 


0 X 

-A = / C? 2_ (MB) 
S X 9^ 




6x 

|r (A+B) - 1-^5 If 


Since 


3f(x-?) 


9x 

-a£(x-5) 
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iiere 



* ^(S-5)^(V 


1 


\ (x-C) 3(,it/2, k) ^ 

® C(x-e)^ + (?„ 


(see 73j[) ^0 of reference [^28^ ) 


Ifence 


5lo o. 

2. = _ [;^(A+B) / - / (A+B) as 

ax ^=-“ 



(x-e) E(V^f ^) 



a 


oo 

n f (A+B-f-D) ’0-(5 j 


(A *6) 
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APPENDIX B 


So deduce expression 


Y/e ha-ve wall integral 

^ - I’ h 

^2 at wall 

= 6 H / (0 + / (^ tn ^r) 3 

(B*l) 

Since a-t c = + 

We haves 

00 i^ <50 

In = eC-/ + f ihr-^0r'> - 

^2 ^ ^c==-y^ 

(B.2) 

ipplyii^ wall boundary condition (3.26) v/e get s 

oo oo ^ ^ 

Iq « 6 C- 

■2 — 

= •-$/ {^0^(5» yyyr) t- 'I' ( 5 ~ ( C ""^yv^ Z] 

+ ?( E. y„) D? ( E =yw) - ( E 3 ^ (®-®) 

(since 0 and are symmetric about c -axis) 
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How integrating by parts the second, e^tpressicn in. above 
integral and applying boundary condition (3t25) we get s 

CO 

C ^ f C ^ C t - y^) •■^4' ( ? — V 3 

(B*4) 
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APPENDIX G 


To oT expression (4 #i) s 




iwSP 


^ f (A-fB+D) t K 

n-^i 


(0.1) 


Integrating each term caa R,H.S. separately! 


n+^i 


/ A 


n-»i 


“ ^p) , y'i dg 

1 1 


wiih a = (Xg«c)^+(^p+^^)^ j b = 4 Ep 

1/2 A Ep r_ 

aixd = CVa) _ = /"— "g- r r - ■ 

1 E^=:X ^ f V TT 1 -t-f V* J 


(0.2) 


(ig» S,)%(Ep 


Taking 3fc. constant over the small interml < E <X.+6. 

J-r 1 1 J 

SO that 0 ( 11 / 2 , 1" ) becomes constant over that interval, 

i 

lisnce (0*2) gives 



(0.3) 
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Similarly / B d? « ^ ^1 ^ ^ * ^ i^ 

^ 2% 


^ i-*i 


In 


V ^ i “ V ^ p ) 


r 


i^^i - g . 


'vT " P - 


^i" ^ i** ^S ^ ^ .* "Sq ) ^+ H ^p ) 


X ^i 




(0*4) 


where i'(^i/2, k- ) and. li('n:/2j k- ) are assumed oonstant over 

^i ^i 

the small interval and 


x .+ s ^ 


p yv ” 

*' / w 

■ y ip^i***^ 




B (%/2 


, te ) i ^ 

^ ^ ^ IlM^wrtP^tinifcii rt -KT^iHwiil fcrTTTrfi^MTi-niifc inU rm^ 


p 


^i 


l + lC ^ * ^***^^^ Jii^ 




1^1 
(0.5) 


Oombining (0*3), (0.4),(a,5) in (1.1), we obtain 


iwSP 


= / Eji Ctan -^ 3. 

2 k * t , ( r ..- 


( r ^. ? p ) 


tan ^ %. 

f - r ^■**1 

^w -^P 


^ hn _3 rllA t^x 


E(il/2. 555 ^) -B(lt/2, l!^^) ^ V’;j|+(^+E^)'2 _ 


In 


2 % 




where % = 

Xg s = Xg - x ^- 


( 0 * 6 ) 
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AEPSKI!XS D 


'io dedaoe eiKpr essioan. (4 .7) : 

^3'*' 

^Xise = 1/ / / 'I'cS 


d dC dv d "> 


(D.l) 


Integrating R.H.S . w.r.t. 5 we get: 


X.4^ . 

i X 




ijSP 


2„ 

•#0 T IL 


r 


1 

" ^ 1 i 
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1 Pd V d P 
3/2 J 

CD,2) 


wh^x-e X^ ~ Xq ^ ^ *•■ ®i 

Xg = Xg - % - 


Further w© hav© 

r + ii^ 

3 3 2^ 

X ^ J s 

<lj o 


Tp + P^“ 


37 ^ 


p d V d P 
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r .+h. 

i 3 „ 

I ZX S 

- O 


dv p d P 


. ^ ^ C o ^ NjJ 


r .'i'h. 

^ 4X ECir/ 2j k) 

"j X^-fCrp +P)^ ll^+Cfp “P)^3 


(D.3) 


(see 291.04 of reference [^28 

xirhere k = / 4* fn 

/ P 

X^+(rp+p)^ 

E(Tr/2, k) Is complete elliptic integral of second kind. 
Further within ttB accuracy of numerical computation vre 


take 


E(ir/2 , k) 




V 2 - 2 

Z +(rp +f) 


constant over in-terval 

J J J J J 

as the point of constancy. We obtain 

^ r.+h. 

J ^ --1 ■ T -- ECtt/Sj k(X)) __X — , 


(X +(rp+rp 


VQj C Cr-p -1)^ 


''^+(rj+ rp)‘ 


1 * 4 * h 

E(Tr/2, kCX) E tan"^ J — L^>± 

X 


•1 r^ -a^- Tt 


(D.4) 



where lc(X) « 



J’inally using (D*^) in (D#2) we obtain 


C (}(2 ) E ( icAi K (%)) 

_ S(Xg) E{it/8, H^z)) ®(Xg)3 


^iiSP 275 
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E 

I. Bvaluatlcaa of Uj^(x, r) for geiieral parabolic body 
of reYQlu.ti.cga s 

The equation for profile of general parabPlio 
body of revolution is given by following two equations, 

The body for which maximum thiolmess lies forward 
of its mid*^oin.t s 

The body for which maximum thicJmess lies aft of its mid 
point 1 

R « 0 H 5^) (B*g) 

max 

The following table gives the data for various G and n 
combinations with position of maximum thiufcness; 



E 

0 

n 

R 


(for max thiclc* (the formula used) 

ness) 


0.3 

i.7i 

6 .03 

(1.1) 

0 .4 

2.36 

3.39 

(E.1) 

0*6 

4.00 

2 #0Q 

(E.l) or (B*2) 

0 «6 

2.56 

S.S9 

(E.2) 

0.7 

1.71 

6.05 

(1.2) 


(see reference ) 
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Ihe value ot (x, r) can be calculated using equations 

(4.2) and (4.3) witb being location of sting ( a 1 

ion no sting). Hcwe-ver integral in equation (4,2) can be 
evaluated analytically for bodies baving maximum thickness 
at mid-^oint* 

The profile is given by s 

H 4 R ( ^ « 5^) 
max ^ 

hence 

S(C)= ISn l4oc ( 5- 5®)® 
s''(C) = 32it E|ax (2 5^'- 3 5 ® + 5 ) 

S"U) = 32ii (6e2_ 6 5+1) 


We put these values in equation (4.2) and 
evaluate the integral of equation (4.2) as follows: 


0 I 2C «• 5 j 


dK « I92n; / 

max 0 


a x‘ 




d5 


■ g^(x- K ) 

15-? 1 


° 1*^ ? I 


Now three cases arise 

i) for S < 0 

I = 19011 1 (1- X - I ) 

ii) for 0 < X < ^ : 

" ” 2 

I = 19011 CSxi @(x+)l)- I + 1 U 
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and iii) for x > 

I SI 192tc ^ ^ -^3 

II . pYaluation of y) for parabolio airfoils • and 

HAQA OOXX airfoils J 


For parabolic airfoil the upper airfoil profile 
is defined by 

(5) = 2 T(5 - e®) 

where is thiclmess ratio of airfoil* 

Putting it in eq^uaticn (3*38) we get 

y) “ f St (1- ) d^ 

It P 0 I (x-O + y J ^ 

^ r (2C-1) In )^4- 

It P >- 


0 


2 «2 


2 f In ) + y 3 


1/ 


dc 


If (to t(5-i)® + C f 3^^" 


" C (5 -^c) In [; (x-c)^ + y^3 -S) 

+ sjf tan'’^ ^ 3 

y 0 

«S , yS 

c> T i iy.1 j i itT . 4- s 

irp ^ ^ (l«x)’^+ y 

2J r tan-^ ^ + ian~^ S 3 

y y 


(E. ^ 
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Similarly for HAOA OOXX airfoil v\re have : 

^ + Bg + Bg 

where B 2 J Bg ana Bg are constants 


henoQ 

» [3 + Bg -f 4" ^ 4* 4Bg 

(E,4) 

First we shall determine the u^C5c, y) cn the surface of 
the Body * Since airfoil is quite thin we can talcs the 
linear velocity perturhation on the surface as linear 
velocity perturbation on the axis y « 0 * So equation 

(3*37) Beocmes! 

V-- = ir { ix -1r- « 

Using (B*S) we get 


Ut(x, O*^) 


1 

^ / J;_ r_™ + Bp+SB,e + 

4- 4Bg (B.5) 


Integrating (B.3) term by term s 

^ dt 

o g(x«5) ^0 (x^ t*^) 


P r 

where t - t 


2 Tr 


U f s 

B- Yx 


0 



6S 


/ d? « 2B. / (-1 + 

o ^ 0 


X •» £ 


) ae 


w 2B^ [[[[ ■*! + S Xn 3] 

1- X 

1 3B^ C ^ 1 « 

f -^dE =3B fC-(S^)-^ ”3^5 

o (x «C) “ 0 S f? 


C- I - 5E + 5® 


and 


4Bo5 


8 


f 

° (x -?) 




de=4Ba / c- (^^ + x5+ 3c^>+ ;:i-.3ae 


s® 


'8 


MM b» 

x:« 5 


ss 4B 


8 


C - i - I _ x" H. # 1« ^3 


JUx 


Henoe (E,3) giYes « 


B, 


o'") “ r? Cl ^ 


l+fx 






3^ X 


2Vx O-tx 

(Bg + 2B^ X + 3Bg S^) 


« 3Bg ( ^ + X ) ^ 4Bg ( f + I + x^) H 

(B*6) 

2be atp-ve solu1:ion for linear perturbation velocity cm the 
plane of airfoil has square root and logrithinic singularities 

at the leading edge (x » o) which gives inaccurate results 
near the' leading edge. It is modified by using following 
formula * 


Ur (5c, 0+) 


1 ul(^» o'") 


(1+ C ^ dY (« )l/2 

^ p dif 


1 


(E.7) 





WALL 


-ij 


3 LOCATION OF GRID POINT IN RECTANGULAR 
ELEMENTS ( axisymmetric flow). 



WALL 


JJMJljLlUJLL 


x x x x x x x x x x x 


X X X X X iX X X X 


X X X X X X 

L-jf: 



X X X X X X X I X X- 


X X X X X X I X 1 X I X 


X X X 

X— L--X- 


(xi.Vj) 

'thi 


4 LOCATION OF GRID POINT IN RECTANGULAR 
ELEMENTS (two-dimensional flow ) . 




L 

C 



U 


( 

I 




i ^ BODY SURFACE ' 

, 1 I -.-i J J 1- ^ L 

0.3 0.2 0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Axial distance , x 


6 DISTRIBUTION OF Cp FOR PARABOLIC ARC OF REVOL 
UTION (M„= 0.90, FR = 10 , free air case). 

o Experimental Present 

(Ref. 29) 











Pressure coefficient, 


“T -T 



Axial distance , x 


FIG. 7 DISTRIBUTION OF Cp FOR PARABOLIC ARC 
OF REVOLUTION (Moo= 0.9, FR = 10, wall inter 
fcrcnce case with porosity parameter = 0. 77 
and wall radius = 1.17) . 

Experimental (Ref. 29) 


O 


Present 








Free air case 



6 

0.77 


10 DISTRIBUTION OF Cp ON PARABOLIC ARC 

OF REVOLUTION (Mco=0.90, FR = 10, wall at 1.17) 


^ ^ ^ « 



1 DISTRIBUTION OF Cp ON PARABOLIC ARC OF REVOL- 
UTION IN FREE AIR AT SUPER-CRITICAL SHOCK 
FREE MACH NUMBERS (FR=10). 




12 DISTRIBUTION OF CpON GENERAL PARABOLIC BODY 
OF REVOLUTION IN FREE AIR (Ma,= 0.9, FR = 10) . 







FIG. 15 DiSTRiBUTION OF Cp ON 6% THICK PARABOLIC AIR 

FOIL AT ZERO INCIDENCE WITH WALL INTERFERENCE 
[M„= 0.825, h/c = 2.0!. 



1 




ii>-ti\<^ (JU.CUbA I’Kb PKKftbHHi;; !iXi>TKJHlJrX(iN UVKH A GI'jNKKAlj 
C [iiJ!)? tJt' KblVLitiUT X UN FkKp, AXH ok In KOHUUis rtAOO 

''SKI. 


AKI, ABOkS AHKJ« 

tiF HulNTfj -JAKK^y on Flibb IiFNGl'H BODY 
Ilf IM'KHV Ai.S TAKh:.-J IH HADIAB OIHtCriUN 

AfT FaXKEXI’ 'lAKFw IN PBKCEWi’ UF BUOY LE.'lbrH 
irtUn MU5, UF iXEKArlONB ALljO<ssf:f> 

K r!ii bung 
r tj c> r 1 N u 
; A ! ! iJN up ;-) i ING 
iiF. Ob BrJl"ib 
IK KHhF AiK CABB 
IK’JAbS, .1 W rKHFB'KKWGE CrtKK 
KAFlO s)F I'HK BUOY 
fjTKBAB (lAu’tl .MiJ 

(MnsiN Nij, nSKI) CAPPENOXX E) 

AKbl CliNft fAs^l'US FUH THE BUOY USBU 
tiUH Hi- NTNO FOMNKli wAUL 
iJBTS’X PAHAOKTKK Ob b^l.li 


JH -.lAlN pRiJCsKA'*! 

1 

t K M , P 

^/l.ri'/B,F,XC!»13 »Kt 3 0) ,HC103 ^I'EXi-OK^KCtsy, 10) ,RBCby) fAHCsU^lOJ 
''s/bf.l/ 1>K#sjX J ,UY^ , Ift.T,fviP 
1 v/BKli/OEbX , IB,BH,1H 
t 'i/BE f S #HK 

t J/KF/K, FK 

I ).lw/I .v/A'KAsJf -H], 10) ,ilU4U; 10) ,KPfRw 

n ^nn/Him/LO ? po , I EW 

a'^io J til.L. I fiO , 1 0 J , Hi. r bu , 1 0 J , IJNbl bO , I U) , U t bo , 1 0) 

A I nSO,.lO) ,EC JOl J *Ft' cion ,0iU40j ,CPL(4U) ,L'PNLt40) 

t r '^bo 

K -lAJ ( lUX , 'TlPb HH NK LV MAX IN XSl KS lb f K H Ifc 0 CJ PU') 
-ppl f. fsH NH IP M '' IN ’^Sl “N Tf pR H IPO Oil ‘O 



m J 


A I’f U»A , iJ'K K'A .^K" ) 

ALO.PT 

f XiajX.ST^Kfi 

I H f i i I jj t* j. i'^ 

V Hi! ^ iKa, CM, P'l 
L , 'V , u; (■•> H I f>; 'f 4 1) ^ H *fli ^ HU P 

)Ai'( JX^ 

lA'i:(/tjX; j, 2,(5X, "i^Hss'^ii^bX, 'LP-', ,12, iX, 'MAX»\ j 

lAl C«T>X ( " t as:' ^11 ^fiX , ''xai’s:^ , Fb,-3, ix, J 

'■'iAi-(iiX ^ ' H.K'', .Cl , /X, ,rb« / ,4X^ 'MS' ,Kb,3] 

.''Al'cax, 'lKas%jit ,YX,'Cus'^b'Cu3#2X, ' PUS ',rb.i 5 
5.S, Mfe'.u} phtt'Uf yili 
ifi,KU,lO PHl.'^fl' V20 

’’A »;( / lux , 'rtiJUX 9fl‘l'H btlNL^*') 

/•A ('(/lux, HKK AIR CAbie<'3 
^■•‘^J‘(/ i(a, "fe'Al.t, 3 fi rp.Hf KHKWCh CAbfcr) 
sy * ijtH 
b V 

b/f- K 
« 4 

- I C iUA ^ 4) 

'I 9 / PUlJAtU^f'3 

U'+ t « 3 /H<ei 

■f Ult 

t j>f(ii<;L,X 

lS,bK»U3 lHS¥^t< ui/ATCj'P^3+tUb,l 

' + X "f P 

'J AS 'at 2 

I '.a J P 
S SSM.iJ 

.SsK'^i 3 Kciwctijfflxcjat ) 

S’! I j -*1 3 tUJ- ?a 

X u U {'. 

f 1 5 sA' u-a 3 +o>'.i(X 

b'f! lu ii 
tear) +0, b^UKjjX 
"■ai 

i( uKUX + Uf.J/X, 

5fji''JjX«'UK 

U JSXC tb I'**! 3 tU« bt CMKbXtOjf t ! 
i +1 3sXi lai ) tux 1 

\ )4U9b'^0>XXtUX^) 



i 


ii- 

1 A *. 

i s 

1-2 , 'JH' 

^ K X 

) « 1 1 +• Ht (.1 ” L J 4 M C 1 J ) 

J e 

J s, u b. , IJ) 11 1 iJ H 3 S ( K *^1 " c K t t'J R « 1 ) ’»' H ( M H " 1 J ) 1 ¥ U , h 

! , I b. e ij ) K ( Is K } s k t ,'J H «* 1 J + 1 U H “ I 3 + H {, fJ K ) 

iKt i)K 3 ,t>.U 3 

M 

XT J 
TsIJid 

“Mt.c 1,1) 

■I!1 t i , j J =^^111 (. I , 1 3 
I “ y* , K 

Wb AK t X !. , K J , Ul,l ^CkL,) 

! ) ^ tT Ij 1 
"IlM^CUKh 
~M1 jU f>J) 

iJUl, Jj 

(t* 

lb 

u 

' .'1 

1st , M A ;t 
.-t , *'iA 

s 1 , t-i N 

I 1 f>’(U A J 

J . AfJJJ, t t ) } CAJsb U^i,h 

'b,i)] CAfib UT ^nTA{Aii^b) 
p)/ )» l-4lb'T2(bb'4 

t PJ+Ai^ib) / (I ."ARU2')'iJl ta,P) ) 

' J sA 1 b 1 '■'1 If 

I , 
t 

-• T , " X 
ST , !lk 

ut, A ,P)-!n f J f lU ) ! (-jj ru lOu 

lb 
lb 

J. f W 



k 


lACC ItiA, I 2f 'TH n fCHAJ iH.i' ) 

l-i ,i'‘A 

U s: 1 , i'i t-f 

, .ns:U ( 1 ,il ) /cma.s 
jUk 

. 

'U n 

!■: 1 f DKi 

k U ^^ULIjU , 3 

f-l. I 'J « fllUhl 1.^1 3«!>RH 

'I (i H, 

-I r 94f LH( JJ ,tis;2^,NK) 

f''’Ar(/22X , 'RSHilOV .StlKr^lU , *'R»'',Ffo. l,iUX, 

J, lOA, ' X ' , / S , ' UL % bX , 

, / A , 'III," , bX, / X, "UU' , •'llwj;,' ) 

, IX U ) ^ I Ui.i. l 1 , J J ^ IJ.^hl I ^ i.i) |i ja= I , ftiH J p i-3 ^ WX ] 

in '34 

nii3u £ / iOA, ’’Cflf' FF, til’ t'HFSb, (Jlfii', !J\) isUUlf WIIHF AUf; *" / 1 4 X , ' X 
X , ''t:i>'>it J KFAiX ' f fix f Ir n'JFftH ' , , '‘Kli'u , FFKT;, '/Kli# tJ J 

I’ 4b, U XC 13 ,C3''{jt U f 3 #Utl , 1 3 3 f 1»1 ,WX3 

K ■ UM- C 1 0 > , f 7 . 4 , bX , f H , 4 , 4 X , (< 4 „ 4 , 7 X , f 4 , 4 J 

<3 l»l,<'JA 
I , S ^ 1 ^ )i f.' 

, J J=fiU , J J 
U J = U]U nH^sJU UJJ 
I « S 1 1". 

J .J I! F 
i jifK 


.’iUHKtiii ! CAbCUFA rKb )Rfc, HOUY HAlUUb Ai^H 1T5 UKH X V AT 1 VfctS 

J'.rslHMJ JjL!CaX1UIV,S 
!^ui) i' I -lb, MjkF £, X , K , ilK , UOK 3 
-Mt 'J/.bF j£ , KS 

lln.-lli''t/M!bS/(:j,H, It'.lj 
~!i . V 
^ H 

■<3 K K 

^ ^\k\\ ) K^.IUUN 

N 

X.UT^ Y } Hh'VWHi^ 
i 



" b- 


u 

uT.o.uuui) (,g TO 111 

C l , -AB J 

L'HVE>»tP"X , )»AB/B 

’IM 


■UwCriUN L'Al,.CtnjATfc;.S DfJKJ, VATlVfc. UF BUpy UKUaS-bECTlUNAl. 
T AX lAli I.OCA'I lUW X 
itpy uSiXj 
SUKFlX,fC,!jK,UDH) 
ii-. VDK 


FUwCl'iLlN CAIjCUIjATEJs SECUi^O DERiVA’lIVE OF BODY CROSS" 
lUNAJ, AHEA A'l X 
T,LC)W U06CX) 

I syKrCX,H,))KrUi)H) 

■•U , ( UH^OK + HN'inlH ) 

IHN 


(;4:4;4;4(sj;»sK;p4(4;3l£4:4i4:)(!!je)(t)(!»4c4r4i4;4s!(£4;»4:)(!s^4:4;;l;f,)l!j(!»^»!^)f;»^Hf3t:))£!l!»!jS4c^si(!)([!!(S)|())C)li)(s:l: 

lb AUXlLlAKl' FUNCTIDm iW C ALCl) B A T 1 nG UHCX) 

CTiiiy f'CTCX^F) 

su,U 

=AB 3 tX-EX 

y IF, I . J T , U , 0 u 0 1 ,) H E X 1 1 K IM 

USl X J-UD£j(E) 

sh'/UiF 

U H H 


^ If ’f ^*** 11 ^ ^******^*!¥ ****** 

i H iii'iCiiUi'J CAliCULAXES UINEAH FEHi'UKBA'XlUN VEBOCIXY OVER 
t .’^UKFACE ax AXlAb LUCAXION X 
JXCh.m UIHXJ 

'1 OH/ b K X I / 1 .) E L X , i B , B B , X S 



tt 


Ci.iM,-1(!''i/SET2/WSflAx, X ,K)S 

UKliSUKliX 

btJM-0,G 

fc. i = U , U 

Uii lU l.=:l , I H 
t-l=:y3 

h,^ = h,1 + U,,5=«^Us::1. 
tisiiJl+UKr, 

hi hhskCt c X , h,i j +rCi cx , b;^ j »4 , u+b c i: ( x , b;3 ; 

hilKs.'iUiVi+SiHp 

L! L J I'! r 1 LI b; 

bLiM = 0b:ii/6,^'hiL’i 

J cju TO luu 

S ijM ss UM t ( X -b, J J / » c b c X ( X , X J +b CT I X , bJ i ) J 

c u w x I w u b: 

hublb'KUh I Y J / ABh ( X- X ) 

(- Atii. SU HF t X , H , IJK , IH)K J 

IF C X . OT . 0 . U . AND. X , I.X. X ) SLEb-hL.KFXDUS (, X J ♦ ALUG ( aH*H*K/ C4,*X*CY-X))) 
UB = bUf1 + hGEb' 

HKTUKN 

end 


XHlh htlHRUiJxlNE CAlliClIi.AXbjy blNFAR FERXUKBAl’lUW VKJjUCI'XY 

ANAX bHUi^ Bum hUKFACF 

5 U H K 0 1 1 r I N K L I N E A H I X , « , U l..i , C' p I j ) 

t L I M .«! U N / a K X .1 / D E Ij X , I, M , B H , j. S 

CUMMlJN/SEl/ Dbj,uYx,UX:^, IhX.wp 

L? bj i ( s D E 1,1 X / 2 • 

Ht< = K13J|‘R=bR 
b; = ”U , »»UF;|j 
hUMsU , 0 
DU 4U!1 J3l,i2 
b.,=:K + DEI.< 

DsX-H.-U.'i^bJElj 

pps I , /.SORT t D + D + HK J 

Nij=i. . /.SMR i’( lD + DEDJ*t'D + DElj3+MH J 

h U 1 1 — ‘ 3 i 1 1' i + i ) a C b< ) ^ ( P P "• D U ) 

t; u .‘j I' J. 1 1 1 i b, 

I b c I.h . lie , JJ GU XU IIH) 
b;sb: + ij , h <' ( Ijb;:i, XDE J 
!) = X"'b>l) , 'S»Db. 

h > I 1 ~h 1 1 f- +. inh t E ) » c 1 . / au R T C D* D+ M H J - .1 , / S R R X ( (, D +U E ) » ♦ 2 + H R J ) 

C 1 • N r .1 M Ll K 

Ui.j = hLiL! 

{:pi)S>'X , ’biji.i 



6 


&u 1^550 , y 

LMJ 10 ISljrUl 

fc-l-yj 

E-isKl+UKl, 

J>lf'iFSKC'i'CX,M Iff C I’CX, II..2J >^4,0+1* etc X ^t:3i 
tUWf li'JUf. 

If C iCi,jiy ,u) Ciu JJ lUU 

bu ’-iJsSU ^ t C Y-t,i 1 / 2 * C f CT C X , )t 1 f f'CT t X , Ei ) J 

CUWfXiVjK 

fcLEb'Sijigj ( YJ/ABSCX"y ) 
eAiUb 3tJHFtX#H^tJK,uyH) 

if ex^BT,0»y«AH[),X,ljT, Y) Sl.fc1t=&ljEKfpUSU)5ftAljUrjCfcSH4‘K’*'H/l4,'^X'^'CY'«X) ) ) 

c, jD 


IHiB bU-BRijUilNE CACiUnpATES ElWEAR PER VU KB AT lUW VEl/OClTY 

AivAy Body bukface 

& U BHOUT i E IjI j'J EAK I X ^ R ? Ob # ePij 1 

tU^i<^yo/SEil/OFbX;lB,BB,iS 

Cur;^tJiv/sei/ D¥.,WiL,in^f IBT,hfP 

Lihl.stJRt, 

hh-hh^r^r 

UPJj 

0 

Oil 400 0^=1 fi2 

haE + OELi 

i}sX«f>0«5¥Ut;0 

pp»l ./SiiRtCU^OfOHJ 

t4M=l , /50R r c tOf Df.bl * C 0 + i5Ey J +BK J 

&uH"5LJMfs)±s(E) ^(PP"*UU) 

.'J.^TiOOE 

lftiX*eOE„Jl 6U ru lUD 
Es=!K+ 0 , b’i'CDEU + uE.) 

OaX^E^O * , 

p/SOHl C C OfUEl^^^fOR) ) 

CU'»J I'lf^UE 
Ut.si>U,^ 

eFL-**2, ’f'Ub 





i'MXS bUHRUUtiNE EVALUATES ThE FIELD IwTEGHAL 
bUBKlJUT Lj'ffc i ftTEGC A i«TUD 
IfifEGfcH $fP 
NEA1.I 

CUrt'^U'M/ J # K tlP J jFtbO, iU J ,RHC^yj » ABC5U, iU) 

COinriuR/bKX/ UE,U^1,DY2, IbTjNF 
CUf'^rtUW/.-aEX.l/OfcliX , IB^BH, iB 
Li li^-'iDfJ/fc.F/K , Ffe, 

Atbi,\OJ#BHaU3,SFilClU3,EClQlJ,rEClUX) 

WAi=L^4-a ■ 

BsSUHTl HD ) 

HKF=:b¥HtP3 

XKCt‘,Fa,l,3 WHPSii^HntS) 

Uij 

bKCJ)=B^HiD ) 

CUPi'jli'JUE 
bUe'^ = U,U 
PU 10 Islff’iXl 
Ps.ijT"UtJj5( 

AFCiS.HE,!) GU TO B 

i f C I . ay , I bT s UK , 1 B Eij , C J.ST + 1 3 3 UFLTbsD 3f i 

XFC-l,EW, CXi>Tti!) 3 l>EjjT®Dy2 

CU.^lTiHUE 

Abl2KlSJ^X{, U+PbS’^DKLjT 
i3'El,EUBlJ GU" TU 2 0 
eai^b^RHt j,«i 3 

HPLstiH,l + BKP 
KkFjshpL>v:RPL 
AH aSiJHTtXbl^XSlfHRP) 

J’CXbH^ASb + HHP) 

AASi-J , U^HRP^bHl 

KRSAK/ C AR=<®XH3 

LAOL iinPULlMs jb#FK) 

HK = AK/ C XRM^XKb) 

CaTjL iRXPUijtKKfEf KKM3 

H b 'I ® BK 3. ■* fl a F 
BHH®2 1 ) “8KX 

(bHK^HM^O 

Alb®BHI^5(EA/XK'tArAOCBMR!S'Xiil/CX5l^XBlFKRi^) 3«EftM/XRbii5A1,’AW{ BHR=S<X3« 
/ UbH-^'XbB + HKrO ? J 
Ab(l“l ? l 

bUWsbUis4i"KFi3*l#1 J^A!Hl*l»lJ 



m 


-J /I J = 

-'LUSSBH(J)+HKP 
'<sSfjHX I X£>H'Aal fRPJ.IJS^HPliUS J 
f\S4*v)*HRP»tiH(,J J / (.XH^XRJ 
Abl, liM J’PUJj|;KK.,K,,hK) 

«?,MysHK! aj«HRH 

Ui ,iJ J = MHt J j/jiR¥EK¥ CATAf^C CRMNS + riHC JJ J/XSIJ-AXANC J/XSI3 ) 

U' C UtQ, i i GU TU 21 

AH ^ 1- 1 , j -] sA U mJ J * A C A « 1 j a ) 

au'^=:,'3U.-i + PU»l 

CU-V 

AA''J) U-,b/J, l4l5y2t>b4^5Uf'i 
K ij. I M H ft 
t 'iU 


iHia aUHRUUTl^E If^TERPOLATEa blNEftHbY h FUf^CTIClW BETWEEN TWO 
PUiftfTh 

auBKOiU’i .! HxpuiiC, AK^FCTj, tf AU) 

U l.iVsKJyalUiV Ecruulj 
CKK-tua,+AK 
«!»c^i^ + l » 0 
UAK“CKK»FbOATt.MJ 

VALsFCTf M,)ttFCTt 5 ^tl)»*rCTCH] J’^Dir 
KLTUHf^^ 


AHl.S b'iHRijUllftig EVALUATES CuteElCiEwTa KUR WAob XHtE(^RAL 
bUbR.iUTl ''(g ^iALu 
i^i'EGEK Sf^ 

KfAb AK#AJl#Mi 

CO‘t^.U'‘i/1AVA^K40,4U,10 J jU1(4U,IUJ,^P,RW 
UiJi'-IOs^i/aEXl/DEt.X , iRj, BH, iS 
■jy '^f -lu^i/Sgl/ Og,OYl^UY2# iST^Np 
CiJi'’'iD4/EFVE,EE 

CU'i.bO,'vj/iXT/K,RfXCbl ) ,K£ 10) jHClO I , NX , NH ^ F C bO # 1 0 J , RB £ $ 0 J , A B 1 5U | 1 0 ) 

Ul^g'^alUM EUOl Jl ,.f'EUUU 

B»5UR1CBBJ 

Brt wsB w 

BKPsSfi^H IP3 

1 ^'tR. EOo 1 JHKPSfj^RBCB) 

BRb-BHW f BHP 



iw y 


bi 3:X t S J -X C 1 ) -f-O , 
iffidaKTtx&x'f'X&i +Hfc<'fj¥aFt»x 
t‘C!5Hy<,ljTipO.MyCUX GU TO i?. 

‘^!MS;bK>!j*®,nKP 

KFJ »/3 si 4 lby 2 &b 4 

^fShO/y I 

is:f^‘» a-U-XSJ =PKI/ 0 !^iO 
lasftL.tiGC-XbJ+YT} 

;oiKAljfJGC 
*U 1 U 
•i-l+i 
j«.L I’siyKtjX 


ti- ( 1 J Gu TU n 

C i 1 # fcHJ, iS‘j' a UH e .1 X • KU* C J. £»T+1 X ) L>Sljf~D Y i 
i)« c u afcw, I ibr+ 2 ,r) uti»T»uyx 

JUU fi 'siOil. 

Sh-J iwx t;>5’*Xtll )+U»S^'D>tbT 
jf i laSUH rtXfoil¥XGJ.l+*BPl.=J^aPb3 
M l=:i,^tXSO/y J 1 

i 4. J » A r A''4 { « XS I X rvi 1 / BH N ) 

AJjX JiKA0U4il“X&Ai + XXl. ) 

CuTl«AOUGt Cla-^KXlJ/tiaX-XlDX 

Rkst, vBKfe'^'BKfV ClXlSj-^XCl) 3 ) 

Cfl.hu i y r poG ( K ft , i*. ^ fc i ; 

OAGIj x.-yTPUiaftft^i' J^^fA j ^ ^ 

1( (tl‘A-,')S(.iU),ttT.«.Va«Olll) ^-lsll.£>»AllBtitl6,0/tS.-K^)) 

A.Ul,s,l/)=UI/l.bJ»Ca.»[Hl-i'IJ-"P'»*CCl.Jl-Cl,I)H2.0/BK»»tti' 

*(AUf 

i twy A i 

K i sift X X 


TTSl'J 1 
COX^CIjII 
A L( 1 ® A ) < i I 
CUN 

oo TO 21 

LON rXNUft 

0 U 20 t“.l#AX 
nsAT X 


OKOTssuBIjX 

LKt IS 5> 'XeJ# A ) CO 'J'O 7 ni,i 

if'tll .KW# J&l tUH* A1 9 tGe ) } Ofe.0T»0y 1 

u CASTA 2A> nfc:OT»OY 2 

{;yNrj"o!K 


X^ll^XCSJ^XtiD+O* 

\ li=!,y<jR l‘( XSlAAtXSiJ 

A* C i f B? P 3 sC C 


/YI j-XSAl/Yll+XSX/YX) 


El) 



1 O'" 


- r ax;->x I 

i 1 

jrvf ri Ok 
rfOK--^ 

JO 


ils bllHRUUTXnsK CALCULATES THK tf^ALb X’^TEGRAL 
jaHUUTl 'XK 1ft ri'^At A1*^L} 

■ilEGbR SfP 

.JM,vtuo/l,^J‘/L»F,ACbaD ,KC10) ,HC10J ,UX,NH,r£b0,10) ,KlH50J ^ABCSO, 10) 
Cuo,v!i) A>^C4y ^40, 10) ,01140 , 10 J 
O -isfO.O 
* ! 10 L " 1 , j'! X 

U4 = KSUH + nUA ,WHJ^A*'JU,S,P) 

UiM r ii'jdK 
X t'i L( ^ b I f ’H 

\0 


ULiiU-nwG AHE the values (JF CUi^PLKTE elliptical iwtegkals or 

Atjj) riKSi' MNpy 
aiAJCK llAtA 

UlAfelLI'jU EUUl) ,FE( lUlJ 
rufi^4U^t/EF/E,FlK. 

OAT A E/i,t>/07 06,Ubf>6B62, , 5bB 948 , 1 , b549b9 , I , 5900 7 i , 

J 1, 3 40 90:^, li54’<^910, 1,518893 , 1 , 5348 i3 , 1 , 530 V 58 , 1 , 5'^b955 , 1 , 522555 
21,518428,1,514284,1,510122, U 5 05 9 4 2 , i , 5 0 I / 4i , I , 49 752 b , 1 , 493 290 
31,489015,1,484/01,1,480456,1,4/6152,1,471817,1,457462,1,463086 
41,45 8 688,1,4542 69,1 , 449 8 2/ , 1 , 4 4b 36 3 , 1 , 4408 >6 , 1 , 4 3 6 366 , 1 , 43 1 8 32 
51,42/2/4,1,42269) ,1,418083, 1,41 3 45 0,1,40 8 791 , I , 404l 05 , 1 , 399392 
64 ,394652, 1 ,3898 83, 1 ,38 5086, 1 , 3802S9 , I , 3 7 5402 , 1 , 3 7 05 1 5 , 1 * 365596 
71*360 64b, Ie3bbb61, 1,350644, 1,345592, 1,340505,1, 335482, 1,330222 
81*325024,1 ,319/88,1,314511,1,309192,1*303632,1,298428,1,2929/9 
91,2y/484,l,28i9 42,l,27b35(J,l,27U/U7,l,265t)]3,l,259'263rl«2534b8 
11, 24 759b, 1,24167 1,1,2356 84,1,229 63 2, 1,223512,1,21 7 321,1,211056 
21*204714,1,198290*1,191781,1,185183,1,178490,1,171697,1,164/98 
31.15/78/, 1*150656, I, 143 396, 1,135998, 1,128 45 1,1, I 2 0/41, 1,11 2 8 56 
41,104? / 5, 1 ,096478 , 1,08/937,1,0 /9 i 2 1 , 1 , 069 986 , 1 , 0 6047 4 , 1 , 0 50502 
51,03994 7 , 1 ,02 8595, 1 .015994, UOOOOOU/, 

IFF/ 1 ,5 70 / 96, 1 ,5 74746,1 *5 /8 / 4U, 1,58 2 780, i ,586868,1 ,591003, 
21*595188, 1*599423, I, 6 0 3710,1, 6 U'8 049*1, 612441, 1, 61688 9, 1,62 1393 
11,625955,1 ,630575, 1 ,635257, 1,640000, 1.644806# 1 ,649678,1 ,654617 

41*659624,1,664/01,1,669850,1,6/5073,1,680373,1*685750,1,691208 



. 7u2i?4# It /y80H7, 1 i.72b77&, t »73iafe*>, 
6U /iH0bb,l*?44i5l,l, /bu/b4, A,7a/^6y,i,7bia9H,i»7?yti4/,l, fTfb'dO, 
/I * /8 4515#r l*7ylbSH>,7 , /g87iy, i»Bybi28,l 1 l»B294feO, 

dl*aj /4yi , 1,8 45894, t.sb^y/ 5,1 ,ab2d4a , i , » / 1 40U, I , HBO i6i , I , W89533 , 
yi,«y8925,l«9U8547,i,9i84iU,l,,92Bb2S,l,9i89UH,l#94958B,l,y&U52l, 
1 ) .'^717 8 3,1^9833/1 , i.9953u3,2,0u7 59a,2,u202 /9,2,0333t>9, 2,U46894# 
22*UbUf?«2,2.0^5383,2.y9l>i/i,2,l9594^,2*1L221i2 , 2 » i 38979 , 2 « 1 5b51 & , 

J2,l/4a7/f2sl9i9/U,2i,214O22,'i.2 45UbS,2,25/2U&,2t28y549,2,i0b232, 
42 « 33 i 499,2,359254, 2 *389015, 2, 4209ii, 2, 45b3iH, 2. 492535, 2 »5iii35, 
5 2.5/8092,2*52/ / /i, 2.083551 , 2 . 7^70 7 i , i * 820752 , 2 , 90 83 37 , i « 0 1 fci 1 1 2 , 
3. IbbS / b,i« 3541 41, 3.595357, 999999, 9/ 

H.' 0 



a^ FKUGKAm LAi.Cl'SjATfcS FHfcfiSiiRt!. MlSTiUBUl'XUK UVEM A I 
1C AJ.ni'i.Ut! iw ifHfcK; AJLR UK XN inJKtlUo *«ALiC WiND 


4Fur VMKiAillibfeS*- 

5S.s0ij.uC culvert? TAREPS qfii AXKeUiL C'MUKU l*fc:>M0iTH 

‘<s?,o&.iJF t^tJA.-*iTh TAt^KN in hA]Cl£,RAh lilRiCi'lU>'t 

i^si-UKC A'^U APT CATK:*rX TAAfe’.S jL^s % Uf AlHFijilj l«fc;piCJ'fH 

AXs,'.|AA j-HUH f<us,ur ITSllHATiUiy^} AULUtsiEU 

s'j'-.UUl^’^Fbij HATiU UF AlKKUlLf 

SFKKii MACh ^ 

tjSU^t'DH tRt.t: AIH CASt 

Sf^UK «ACIj lw'i’KHFS?;HEWCfc. CA&E 
ftSriAC? HMiiHT U¥ I'Hfc: isfXE^U I'UWNKLi 
CWCUrtCtal i'5f CARAHE.'A"fc;R Ofe* 


I%ll> if' TMK J*^Ali5i CHiJGKAW 
t.Au ‘’i 

icrh* f.'f s#P 

-J .‘'1 4 ' / 

lUi'i ‘Ht.j/.i2/Ux (40#1U) , rflAUj-IO, iuj 

L'iP^('jn,,/y3/X(40? # YClU) # H 1 1 U ) » WP ^ , UEUK , ¥# 

;Un'lu./!f4/F 140^10) , A(4U , 10} , (40 j ^ X (40) 

, O "'. ^ ' i ''< / S' i( / X 

:oM 5i}'!/Ul]/UcL.(4U,XO} ^UNCCfU, IM) 

Ui,4U,Xf>}, UliC40#10} ,AAlXu) 

i!,KK JhS.>.0<ti 

T«Pa ysu 

a -.jy^AT C J ox , WB r f CP# « AX, T# # X U ^ ) 

Aeof.Pi' NB,NK,CPj.MAX,T,M, XC 
.11* f I '..s^jiijO) lypfc. obi 

!• j). ..-41 (IDA, *T¥Pfc lf«^,WP'} 

li; { n'j.fXKeU? ^ACCe.PT 
P<M'J a30,Nb,wR,0F,>MAX 
Fsa-c-s BXO, Iti,'i',M 
IKXLC.J^F.O} PKliiJT 8i0,xw,wp 

ruit -Af {/6A,’*flfS«'*',l2,fcX,’'RH«''#iy,t5A( •'lfPs‘^,l2,3X,''MAX 
KuK.''i'iC^*X,^XL»'’,lJ ,yX, 'ts'',fb, 3,4ft. 

S' UR (A, (OX, ^XR»’',PS#3,3X,'HFs^,i‘S5, Jj' 

It- £ U,,FU.U) PKiW'f 910 
U C PKlf'*T 920 

Fuh 1 /- ! (/iux, 'FHfcE MR CA.SF' ) 



^ »» 


i / X V A f 


AiLi* I i-Zi-’r. Cr; 


isi 

A 

r/ 2 s 

i a 'V 

,LX™s a /rijUATl. wH) 

■ ( 1 , ) 

.1 jrt u, 5’“1 ''JjOATC«FJ J'FJJfc-ijX 
f 3 <.' 

IJ 3 SA ( JL«J )+yE;i»X 
Jii J X f'iUS, 

f ( J L,‘yt!., 'U GlJ lU il 
y I, I J-f -lAX 
i/}^ nhA 

► H X 3 s i ( I 3 

UU /O jssJfftJK 
M Ci3 3-MCd*” i 

UM’X’vU^^ 
iL' I’C lA- 

‘ H^U t l^k^t } 

-( I JsruH/2. ' 

( V X“Hf 3 > 

K2 J»L3/3 

i U 

3U 2t JssJ,f"3K 
1tU3=HC.(»-lX 

UJX»?f!j«XX+HXJ”XX+HC!j) 

ilL=^#KX3^HU 

L'Ui'*) 1= i«Ut. 
tXU 11 
S 1 “ ^ t J X 

Jf-iFUl 1 )=¥UAiX 
yj»u«uouui?i 
l i ? 1 3«CUtMfe-*<'UtiXXl/ X'JX 

f t J. , 1 J-lJ-iUf 1 U 

Ulti, fXX-USUrl) 

UU iO a=2^«H 
S s > c 1 >^n 

U“ UW&^33LUlr¥4X 

f i i , J X'siiB UfJ)^ViiLl,J3 



^ ’ 


t c U !btl , i) 

J llf L f' 

1 rr i i* 

( nE(>„ 0 

i b“U , (' 

> ‘iO irsi,5*Ax 
ff ?s?3 ,|.jX 
n 42 P"i^fV'Y 
C<M,b i ‘.ITKGCAIrJTfO 

IP’C r!\,KO, u AryfJ* IL.MI=:,01 CAIjb ^At.L 
TFCTI .fs'E.O) CALI, W^hlM'^U} 

ACPl^AIfjTG+AlWL 
• 4ijy=;A CS , P J 

SIC.‘^^^r')rCHBCS,PH'AINTG“Awy2^F(S,P3 +Ai.WL)/C J ,«AHU2’^lU CS,P) 3 
’fJf'iTT --vsiJK 

>n bi f^lfhx 

!>n JSS!,HY 

th (aP-,S(UCI,U)®U3 C, 1,^133, fiT,ERHEi|?) GH TO 100 
3ur:Ti *nK 
ry . C3 

pKf-n 9 2, IT 

F'-iRf- Aft 10X,I2, 'TH ITERATItP^^? 

PKI-^T 93, CYC J) fGasI ,NR3 

FHR iATC/2iX, 'Y^'',Ff^,3,10X,'VK',K5,3,lOX,'¥s'',r6,3,iyX, 

J '¥~%Ffj*3,10X, 'Yss%Fa,,3, lOX, 'Ys^,Fb.3/14X, ^X%7X,'UL% 
?bX, ,7X, 'nL%6X, 'USSL' ,7X, 'MI,' ,foX, , ?X , 'UL' , 

3<sX, ''U''JL‘',7X, 'UL',?iX, 'U*yL',7X, 'ML'#6X, 'UNL') 
r,f) b? 1 = 1, MX 
bf! b.i u = l,My 
i M,CI,M)=UCT ,.n/C£)MK 
. ( L( f , JJ/Cn^‘S 

cnvT r.iijc: 

PKfiT 9l,X(T3,( MLLC T , J5 ,MNJ.t T , J) ,J=i,ft|R3 
rn‘->Ti iJs'E 
CALL UOTpiJT 
b i E)P 

pf) =>2 J=1,MX 
DO fjt 0 = 1, •'^Y 
I ICI , J)=CU1, J) 
t ( r ,u)“:i Cl r , 01 

Cfi-Ml nt: 

CL*'' f r ' UK 
CL jti- up: 

hTQP 



H- ,u..iin ? .F rAr..r.}!.Af r.E t-lKl-t) T«TL{^RAL 
i'iiiiwR AT^^TfO 

•Vv.,s fi,V 

>' y 'l/S,P,l-4K,K‘H 

''i'^P'*‘/0J/Xt40> ,yClC) ,H( UU Yii 

H! ■/n4/KC4f|,]l03,AC40,10) ,¥BRl/C40 3 ^ VB^M( 40) 

;!■? 

:C,.- 

/.A *141^02654 
'ff .'')-XC3 )+0,5>f;nE!jX 
) fi s, s: ] , A( X 

'/ 3 e* !.'* f ^ X 

X f B 

I ^'J»/i,=S=HCn )^B 

( ' J ^ ft 

~A rA*‘4fYl/X'^)«ATAMCyi/X) ) "ATAN ( ¥ ^/X 2 ) 4^ AT AN ( Y 3/X I ) 
s-AVArj{Y2/X2) + ATAN(Y2/Xl ) + ATAN ( ¥4/ X2 ) -ATAN CY4/XU 
, 1 )~-r=S=CAll + A12) 

]y .)S2,rJY 
‘ / 

’Y4t?s ^HC 

'-“AViU4CyV/X2)4-ATANCY2/Xl)+ATANC¥4/X2}-ATAN(Y4/Xi) 

Jr«-c*CAil+AJ2) 

l- I j H ! _ 

= x 2 
S' 1 ' ’ ' ! h 

f7*“ SSI's) 

ins 'I 


T.s .SOHRUnTlNK KVM.nhTKB CtlKPFICIKNTS FUR «'At.L INTKGRATj 

mFOi<Tf'*lF WAl.Tj 

;TKf;r-;R 

j.. n=fV'i?/in (40# 1 03 ,s*.'t40,40#io) 

Cfjfnj!f^,VO'3/X (40) #¥CiO) ,H(10) # WP^BrOFLX# YW 
inr,i/ii4/F (40# tO) #AC40, 10) # YBPLC40) , ¥6^fii(40 3 
Y s; V* H 

i 41 b 026^4 
1=0'-^ ( Y ( P)«YiO 



/ 1 ^ f 1 

Xl^Xt 

l , rjt;t tXXl 4 yY 13 > r ( xxi + yY : i )) 

TAfrcxl / Yl)«ATAMCXl/y2) 

0 3 .- S 1 , fa 
l-PKI.X 

iLiiKC CXXi!4Yyi 3*!'CXX2tyy233 

\Taiux2/ya 3«A’rANf x?/y2) 

Teu, b=i(’^fp/8#(A2“'Al 3«Ca2"Bl ) 

, S , P)sCWfDI 

<2 

32 

YZ 

ri'-ifiK 

JH '- 


5 r->UhRq[!TlNK CAfjCULATK3 wAbb 1, ^STKG'HAb 
.s^O(JTr\‘E .f NTififACAiy^l,) 
n,P 

?'iO!^/nt/a,p,?U/NR 

' O '^’ M DM / n 2 /Ij 1 C 4 0 , i 0 ) , W c 4 0 r 4 0 , 3 CO 

S'JR 

to I = 1|,KX 

nssiHij '' 4SJ i { .t OJ Y J C X , B , P ) 

in ’ O-oK 

TMRO 

li 


Jfs S.IrtRfiOTl i^T CfibPULATKS PPKSSUK^: DlSTKIBUTIOM 
ifC^‘’nft tvfv nUTPOT 

^! n !/ nii / iiLi .(40 ou ) , n ?^ bC 40 j 103 

oni/Xt40) ,yci0) OKIO) #E«P,B,UEr,K,YW 
HUJ 2 ]=AJ ^ ( tt , 4 S 2 ^C ! »» Q 23 ) 

F Jlt“3AT(tUX,^MACH 3, lOK, 'PARABOLIC AIRFOIL") 

FUfir'ATCSX, 'X CP LI« CP PRCBRNT') 

npMATCbX # F 7 * J , SFll .33 





PRI^T 9^,*^ 

P«lNf 

A'iK(G*-i , 1 *S^/ 2 , 

KXs«/tfi»l , ) 

10 Islisf^X 

Xisxc.1 3 

DiXsDYlXX) 

QbSid^i-illhCX f 1 ))Y^* 2 ^ 0 YK’^DYX 
QH^Cl .“‘■U il.t 1,3 )}4E#2+D¥X«;OyX 

cpr, jtCPcaL) 

CP-j^CPC'iHi 

PRJ.JT 94,XI,CPt,iCP^J 

cau PTotii^ 

RVT’JR^^^ 

Ef-iD 


;ft f J«r 4 #: y 4 ¥ y 4 * « sS! X: Jf! # 36« f ® ^ ^ ^ y y ^ sf; y 4 ^ ^ Jj( -f: 4 Xe ^ y ♦ 4: f 2^ 5^ * 4 * !j: * *: * 4 i}! t 

THE FDL-L0i^If*<3 PACKAGE BE REPGACEf) BY APPRUPIATE FUNCTlflM.B 

IF THE BODY PHOFTIjR IS CHAMGEO 


^ !jt « f $ 4: s^ s)E ^ ^ JSC ^ )js ^ * Jf: J^s $ ^ ^ 1^ « t It' * f * Xs ^ ^ S ^ y ^ )|= y =Jf » J}! ♦ 3(c jjc ♦ ♦ 4: ♦ $ if: 4 4 

the rr»LOO^lNG PACKAGE CALCUOAtE.S RlfRCriD^S FOR BODY PROFILE 
ITS DERIVATIVE AVO LTMEAR PERTURBATION VELOCITY FOR & PARABOLIC 
AXRFOIL. 


This ruocTirjfi calculates airfoil prdftlk 

ro.ycTto-f YYCX) 

CUH’^fPt/TI/T 

YYsOoU 

1 F C K « LT , 0 » U » OR « X * GT » 1 <i 0 3 HFl UH U 

YYs 2 «O 4ri«X^{l»0«X3 

RETURN 

E-^jD 


This kil:^ctidn calculates the oerivaiive of 
FUi^CTlON l)YCX3 
COO.ynH/rv/T ' 

DYSO , 0 

IFCX, L'UA.O, UR, X, 01.1,03 RSTURS! 


AIRFOIL PROFILE 
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ID 


a'i I'UwcTilHi’ CALCULATES LINEAR PERLUKBATtON VELOCITY 

'r^CTTO^^ ULCXjY) 

CnN«nr!/o.VXTC401 ^Vf c lO) ,HC io ) ^B.DELX, Y¥ 

iu *0^’/ rx/T 

- 3, 14lf5D2hS4 

4K:2«Uf‘lVPl/S»C CO,S»X)+ALUGC CX^fX+y^^Y)/ c C S. * 0 

l-fY^O-^.O^^ysjtCATANex/Y^ + ArANC Cl ,0-X)/y)n 
CX ,GT.O s , ANU»X,LT.i» ) tJLsUL+ ( I , C 3 » /PI/PI 
? , , Y ( X« « f? 3 dialog [ X/ c 1 «*X ) ) ) ^ « /PI /Pi’f' C ALDGC X/ C 1 ,«K 
) ^=4^ /» C 1 . *4 , sp C X*- , S 3 ^=*=23 ) 3 
ruHSf 

u 





